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SUMMARY 


NASTRAN contains two techniques to solve the differential stiffness prob- 
lems. One is incorporated in a new static analysis rigid format and the other 
is contained in a new normal modes analysis rigid format. The purpose of this 
paper is to compare the two techniques relative to computational accuracy and 
time of execution on Level 16.0. 


INTRODUCTION 


Through Level 15.5, the Static Analysis with Differential Stiffness 
(Rigid Format 4) capability was based on a one-step process (Reference 1). 

This process was a static solution to obtain the linear stiffness matrix and 
internal element forces followed by an element differential stiffness comput- 
ation. This procedure was based on the assumption that the internal force is 
a linear multiple of the applied load and that the applied load remains fixed 
in magnitude and direction, moving with its point of application. The user 
provided differential stiffness linear load factors. An iterative technique 
was introduced (Reference 2) and is now fully described in Reference 3. 

The new approach to solve the differential stiffness problem is begun 
with the iteration of the displacements to compute the differential stiffness 

matrix from 

[K+ k"* (u,)] = {P} (1) 


where u. and u^.^-j are the set of displacements at two successive iterations, 

K is a stiffness matrix, and P is a load vector. Rearranging terms, [K (u^.)] 
is removed from the left hand side and is replaced with the term [K (u^)] to 
give 

[K + (Ug)] {u.^^} = {P} + (Ug) - (u.)] {u.} (2) 

or 

(K + (Ug)] {u.^^> = {P} + 


149 


estimate initially equal to the linear elastic solution. With 
internal loads may change due to differential stiffness 
effects so that the solution is not linearly related to the applied load. 

Thus equation (3) treats the change in differential stiffness as a load correc- 


Pi^ovided to control the iterative process. The first 
BETAD, limits the number of load corrections before adjusting the differential 
stiffness. The second, NT, limits the cummulative number of iterations. Thus 
load correction iterations can be performed up to the limit BETAD, at which 
time the differential stiffness is adjusted, and then more load correction 
Iterations are performed and an adjustment is made to a new differential stiff- 
ness until NT is exhausted. The third, EPSI0, is a convergence criteria which 
terminates the process when successive iterations of the differential stiff- 
ness are sufficiently small. Convergence occurs when e. < EPSI0 where 


» 



T 





(4) 


The user^either relies on the default values of BETAD=4, NT=10, and 
EPSI0=l.OxlO"^ or prescribes values through a PARAM bulk data card. 

Figure 1 is a simplified flow diagram of the procedure. The requirements 
Of the rigid format are that two subcases be used to define the static output 
requests and the differential stiffness requests. Loads and constraints are 
defined above the subcase level and plot requests are last in the Case Control 


A new normal Normal Modes Analysis with Differential Stiffness (Rigid 
Forirat 13) capability was described (Reference 2) which combines static, 
differential stiffness, and normal modes analyses. 

Presently, this technique is based on the original differential stiffness 
approach (Reference 1), but is limited to one loop (or load factor) through 
the Rigid Format. ^ 

Figure 2 is a simplified flow diagram of this process. The rigid format 
IS utilized via three subcases. The first pertains to the static analysis 
Where the load is defined, the second prescribes one load factor for differential 
stiffness, and the third contains a method for a real eigenvalue analysis. 
Individual output requests can be made at the subcase level and plot requests 
are last in the Case Control deck. 
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TEST CASE 


The test case used is the standard NASTRAN Demonstration Problem for 
Rigid Format 4. The structure is a hanging cable acted upon by its own weight, 
which is an equilibrium position, assumes the shape of a catenary. The original 
shape of the cable is circular. The final shape of the cable is readily pre- 
dictable from equations developed in Reference 5. 

The coordinates of a point (x,y) on the catenary are defined by 

X = S sinh-1 (f ) , 

(l + i'ff) -l] . (6) 


and 




where H is the tension at the bottom of the catenary, w is the weight per unit 
length of the cable, and s is the distance along the curve. 

The original demonstration problem is in English units and are converted 
to the Newton-meter system for this discussion. The input data decks are 
shown in Tables 1 through 5. Notice the alter for Rigid Format 13 necessary 
to allow multiple load coefficients and plots for the original differential 
stiffness technique. 


RESULTS 

Tables 6 and 7 show the deflection results obtained by the two techniques 
compared to the theoretical expectations. The results computed by Rigid Format 
4 were obtained in four iterations when the convergence criteria changed from 

4.5x10”^ to 4.3x10'®. On the CDC 6600, Functional Module, DSCHK, which performs 
differential stiffness computations in Rigid Format 4, spent about 2 cpu seconds. 
Functional Module DSMGl (used four times) consumed approximately 5 cpu seconds. 

The results shown for Rigid Format 13 were those obtained after the first load co- 
efficient. (Successive coefficients produced deteriorating answers.) Using Rigid 
Format 13, Functional Module, DSMGl, used 1.25 cpu seconds and Functional Module, 
DSMG2, used about one-third of a second. 

Figure 3 shows the graphical results of this test. 

CONCLUSIONS 


For a simple structural element case that can be readily verified, there is no 
appreciable difference in the results computed. In fact, Rigid Format 13, when 
using only one differential stiffness coefficient, actually computes the nonlinear 
solution faster than the Rigid Format 4 counterpart. Thus, the two differential 
stiffness techniques available in NASTRAN can be utilized equally well depending 
upon the user's preference. 
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Table 1. Executive Control Deck for Rigid Format 4. 


ID 0IFFSTIF*RF4 

APR DISPLACEMENT 

SOL ^*0 

T I vtP 10 

CEND 


Table 2. Executive Control Deck for Rigid Format 13. 


ID DIFFSTIF»RF13 

$ REQUIRED ALTER TO CORRECT AN ERROR IN THE RIGID FORMAT 
ALTER J> 

DARAM //C«N»ADD/V»N»N0MG6/C»N» 1/C»N*0 % 

$ ALTERS TO CHANGE THE RIGID FORMAT FOR MULTIPLE D, S. FACTORS 

ALTER 13R i 

JUMP OSLOOP S 

LAHEL OSLOOP 'i 

ALTER 1S5 S 

COND LBL80*REPEAT0 S 

REPT OSLOOPflO 1> 

PARAM //C»N»NOT/V»N*TEST/V*N*REREaT $ 

LABEL LHL80 T 

alter 158» 17a $ 

ALTER 176» 176 S 

»LOT plTPAR.6PSETS,ELSETS»CASECC»BGPDT,EQEXIN»SIL»»PUBGV1»GPECT , 

OESB 1 /PLOTX3/V ♦N»NSIL/V»N*LUSF T / V » N * JUMPPLOT/ V » N *PLTFLG/ 
V*N»PPILE 
alter 187* 188 i 

ENOALTER S 

APP DISPLACEMENT 

SOL 13*0 

time 10 

CEnD 


153 


Table 3. Case Control Deck for Rigid Format 4. 


title = DIFFERENTIAL STIFFNESS ANALYSIS FOR A HANGING CABLE 
subtitle = RIGID FORMAT 4 SOLUTION 

LABEL = INITIAL SHAPE IS A CIRCLE* FINAL SHAPE IS A CATENARY 
LOAD = 32 
SPC = 2 

displacement = ALL 
SPCFORCE = ALL 
STPfTSS = ALL 
FORCE = ALL 
OLOAD = ALL 
SUBCASE 1 

LABEL = LINEAR SOLUTION 
SUBCASE 2 

LABEL = NONLINEAR SOLUTION 
BEGIN BULK 


Table 4. Case Control Deck for Rigid Format 13. 


title = DIFFERENTIAL STIFFNESS ANALYSIS FOR A HANGING CABLE 
SUBTITLE = RIGID FORMAT 13 SOLUTION 

LABEL = INITIAL SHAPE IS A CIRCLE* FINAL SHAPE IS A CATENARY 
SPC = 2 

DISPLACEMENT = ALL 
SPCFORCE = ALL 
STRESS = ALL 
FORCE = ALL 
SUBCASE 1 

LABEL = LINEAR SOLUTION 
LOAD = 32 
OLOAD = ALL 
SUBCASE 2 

LABEL = NONLINEAR SOLUTION 
DSCOEFFICIENT = 50 
BEGIN BULK 
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Table 5. Bulk Data Deck 


$ GEOMETRY IN CYLINDRICAL SYSTEM MEASURED IN METERS 
$ 

0.0 


COOD2C 

10 

0 

0.0 

0.0 0.0 

♦COORD 

1.0 

0,0 

0.0 


GRID 

10 


3.048 

0.0 

GRID 

11 


3.048 

10.0 

GRID 

12 


3.048 

20.0 

GRID 

13 


3.048 

30.0 

GRID 

14 


3.048 

40.0 

GRID 

15 


3.048 

50.0 

GRID 

16 


3.048 

60.0 

GRID 

17 


3.048 

7C.0 

GRID 

18 


3.048 

80.0 

GRID 

19 


3.048 

90.0 

S 





S CONNECTIONS VIA BARS 



$ 





BAROR 




-1.2 

CBAR 

10 

10 

10 

11 

CBAR 

11 

10 

11 

12 

CBAR 

12 

10 

12 

13 

CBAR 

13 

10 

13 

14 

CBAR 

14 

10 

14 

15 

CBAR 

15 

10 

15 

16 

CBAR 

16 

10 

16 

17 

CBAR 

17 

10 

17 

18 

CBAR 

18 

10 

18 

19 

S 





S GRAVITY LOAD 

FOR NEWTON-METER 

SYSTEM 

$ 





GRAV 

32 

0 

9.8 

0.0 1.0 

$ 





S CONSTRAINTS 




GROSET 





SPC 

2 

10 

12 

0.0 19 

S 





$ MATERIAL AND 

PROPERTY 

DEFINITIONS 

$ 





MATl 

1 

2.63^7 


.3 1.78-2 

PBAR 

10 

1 

9.29-3 

6.87-6 6.87-6 


I.O 


0*0 


0 

1 


$ DIFFERENTIAL STIFFNESS PARAMETERS FOR RF4 
S 

BETAD 8 
EPSIO 1.0-5 
NT 18 


PARAM 
PARAM 
PARAM 
S 


$ DIFFERENTIAL STIFFNESS COEFFICIENTS FOR RF13 
$ 

OSFACT so 1.0 2.0 3.0 A.O 

$ 

ENODATA 


5.0 


0.0 


0.0 


345 

0.0 


1.0 


1 


♦COORD 

TOP 


BOTTOM 


TOP 


BOTTOM 


RF4 

RF4 

RF4 


RF13 
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Table 6. Horizontal Deflections 



4 


Table 7. Vertical Deflections 


Grid Point 


11 ( 10 °) 
13 (30°) 
15 (50°) 
17 (70°) 
19 (90°) 


4.25 

3.19 

2.12 

1.06 

0.0 


Theory 


-0.00341 

-0.00696 

0.000914 

0.1737 

0.2846 


Uy - Vertical 


NASTRAN RF4 


-0.01245 

-0.03865 

0.04478 

0.2430 

0.3707 


NASTRAN RF13 


-0.01241 

-0.03432 

0.09006 

0.2418 

0.3679 






/ NT \ 
Iteration 
.Limit / 


Input Processor,, 
Geometry Processor, 
Undeformed Plots 


Differential 
Sti ffness 
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Constraints 
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Figure 1. Simplified flow diagram of Rigid Format 4 procedur 











Figure 2. Simplified flow diagram of Rigid Format 13 procedure. 
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Figure 3. Hanging Cable Test. 




